Graviton production in the scaling of a long-cosmic-string network 
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In a previous paper [1] we considered the possibility that (within the early-radiation epoch) there 
has been (also) a short period of a significant presence of cosmic strings. During this radiation- 
plus-stnngs stage the Universe matter-energy content can be modelled by a two-component fluid, 
consisting of radiation (dominant) and a cosmic-string fluid (subdominant). It was found that, 
during this stage, the cosmological gravitational waves (CGWs) - that had been produced in an 
earlier (inflationary) epoch - with comoving wave-numbers below a critical value (which depends 
on the physics of the cosmic-string network) were filtered, leading to a distorsion in the expected 
(scale-invariant) CGW power spectrum. In any case, the cosmological evolution gradually results in 
the scaling of any long-cosmic-string network and, hence, after a short time-interval, the Universe 
enters into the late-radiation era. However, along the transition from an early-radiation epoch to the 
late-radiation era through the radiation-plus-strings stage, the time-dependence of the cosmological 
scale factor is modified, something that leads to a discontinuous change of the corresponding scalar 
curvature, which, in turn, triggers the quantum-mechanical creation of gravitons. In this paper we 
discuss several aspects of such a process, and, in particular, the observational consequences on the 
expected gravitational-wave (GW) power spectrum. 

PACS numbers: 04.30.-w ll.25.-w 98.80.Cq 



I. INTRODUCTION 

CGWs represent small-scale perturbations to the Uni- 
verse metric tensor [2] . Since gravity is the weakest of the 
four known forces, these metric corrections decouple from 
the rest of the Universe at very early times, presumably 
at the Planck epoch [3] . Their subsequent propagation is 
governed by the space-time curvature [4], encapsulating 
in the field equations the inherent coupling between relic 
GWs and the Universe matter content; the latter being 
responsible for the background gravitational field [5] . 

In this context, we explore the creation of primordial 
GWs during the scaling of a long-cosmic-string network 
(see, e.g., [G]), i.e., in the transition from a cosmological 
model in which the matter content can be modelled by a 
two-component fluid - consisting of radiation (dominant) 
and cosmic strings (subdominant) - to a pure radiation- 
dominated Universe. 

Cosmic strings are one-dimensional objects that can be 
formed as linear defects at a symmetry-breaking phase 
transition [7], [s]. If they exist, they may help us 
to explain some of the large-scale structures seen in 
the Universe today, such as the gravitational lenses [9]. 
They may also serve as seeds for density perturba- 
tions [10], [11], as well as potential sources of relic gravi- 
tational radiation [12], [13]. 

The presence of cosmic strings in a radiation model, 
is responsible for the constancy of the effective potential 
which drives the temporal evolution of a CGW, leading 
to a critical value (kc) of the comoving wave-number, 
which discriminates the metric perturbations into oscil- 
lating (fc > kc) and non-oscillating (fe < fee) modes [1]. 
As a consequence, the propagation of CGWs through 



a radiation-plus-strings stage would leave imprints on 
their power spectrum. However, such a stage docs not 
last very long, since, gradually, the production of loops 
smaller than the horizon results in the scaling of the long- 
cosmic-string network. According to this process, the lin- 
ear defects form a self-similar configuration, the density 
of which, eventually, behaves as and the Universe 
(re)enters in the pure-radiation era (see, e.g., [7]). 

In the present article, we study the quantum- 
mechanical production of gravitons in the transition of 
the Universe from an early-radiation epoch to the late- 
radiation era through the radiation-plus-strings stage. 
Several theoretical (and observational) consequences are 
discussed. 

The paper is organized as follows: In Sec. II, we 
summarize the theory of CGWs propagating in curved 
space-time, including the implications arising from the 
constancy of the effective potential. In Sec. Ill, we 
take advantage of the scaling of a long-cosmic-string net- 
work to arrive at the late-radiation era. In Sec. IV, we 
demonstrate that, due to the modification in the time- 
dependence of the cosmological scale factor along this 
transition, the scalar curvature changes discontinuously; 
hence, gravitons are created, and, in Sec. V, we explore 
the corresponding power spectrum. The high-frequency 
part (fc > fee) of this spectrum acquires a character- 
istic profile, resulting in a periodic function of the fre- 
quency. Finally, in Sec. VI, we evaluate the (integrated) 
energy density of the CGWs created in the scaling of 
the long-cosmic-string network and we compare it with 
the corresponding quantity predicted by inflation (see, 
e.g., [12]). It is worth noting that, the energy density 
of the CGWs created in the transition of the Universe 
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from an early-radiation epoch to the late-radiation era 
through a radiation-plus-strings stage, not only is added 
to the corresponding inflationary quantity. In fact, at 
times 10~^^ sec < t < 6540 sec, i.e., within the first 
(and, certainly, most interesting) 109 min of the Uni- 
versal evolution, it dominates over every other form of 
GW-energy density. 

As regards the space-time geometry itself, we confine 
ourselves to the spatially flat Friedmann - Robertson - 
Walker (FRW) model. This model appears to interpret 
adequately both the observational data related to the 
known thermal history of the Universe and the theoreti- 
cal approach to cosmic-string configurations [7]. There- 
fore, it seems to be the most appropriate candidate for 
the curved background needed for this study. 



II. CGWS IN A FRIEDMANN MODEL 

The quantum-gravitational creation of gravitons in an 
expanding FRW Universe was first demonstrated by Gr- 
ishchuk [14]. He showed that, in the linear approxi- 
mation, the behavior of a CGW propagating in curved 
space-time is identical to that of a massless, minimally- 
coupled scalar field. In particular, each of the two po- 
larization states of the metric perturbation satisfy the 
Klein-Gordon equation 



Kj ;a = 0, 



(1) 



where Greek indices refer to the four-dimensional space- 
time, Latin indices refer to the corresponding three- 
dimensional spatial slices and the semicolon denotes co- 
variant derivative. The quantization of primordial GWs 
and that of the minimally-coupled, massless scalar fields 
also proceeds along identical lines, as it was demonstrated 
by Ford and Parker [15], [16]. 

A weak CGW (|/iy | <C 1) propagating in a spatially- 
flat FRW cosmological model is defined as [12] 



ds' 



c^dt^ - i?2(i) {S,j + h,j)dx'dx^, 



(2) 



where c is the velocity of light, 6ij is the Kronecker sym- 
bol and the dimensionless scale factor R{t) is a solution 
to the Friedmann equations with matter content in the 
form of a perfect fluid. In terms of the conformal time, 



dt 

Wr 



(3) 



the solution to the Klein-Gordon equation (1) is a linear 
superposition of plane-wave modes 



/i.y(T,a;'') = (/jkir) £ij e"' 



(4) 



where (f)k{T) is a (complex) function of time and Sij is 
the polarization tensor, depending only on the direction 
of the comoving wave- vector fc^. For a fixed wave- number 
fc^ = X) ^ri the time-dependent part of the corresponding 



GW mode satisfies the second-order differential equation 
(see, e.g., [17]) 



0, 



(5) 



where a prime denotes differentiation with respect to r. 
Under the further decomposition 



hkir) = Mr) i?(T), 



(6) 



Eq. (5), which governs the temporal evolution of a CGW 
in a FRW Universe, is written in the form 



K. 



0. 



(7) 



Eq. (7) can be treated as the Schrodinger equation, for a 
particle moving in the presence of the effective potential 



(8) 



To solve this equation, one needs an evolution formula 
for the cosmological model under consideration. 

In terms of the conformal time, the spatially-fiat FRW 
model is a solution to the Friedmann equation 



R!_ 

i?2 



SttG 



P{r) 



(9) 



(where G is Newton's universal constant of gravitation), 
with matter content in the form of a perfect fiuid source, 
T^u = diag{pc^ , —p, —p, —p), which obeys the conserva- 
tion law 



.i?' 



P 



1 



3^ [P+-P]-0 



and the equation of state 



P 



1 pc^ 



(10) 



(11) 



where p{t) and p{t) represent the mass density and the 
pressure, respectively. 

The linear equation of state (11) covers most of the 
matter components considered to drive the evolution of 
the Universe, such as quantum vacuum (m = 0), a net- 
work of domain walls (m = 1), a gas of cosmic strings 
(to ~ 2), dust (to ~ 3), radiation (m = 4) and Zel'dovich 
ultra- stiff matter (m = 6). For each component, the con- 
tinuity equation (10) yields 



(12) 



where Mm is an integration constant. Provided that the 
various components do not interact with each other, a 
mixture of them obeys [18] 



P ~ 2^ Jim ' 



(13) 



where, now, Eq. (10) holds for every matter constituent 
separately. 
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A. Constancy of the effective potential 

There is a case of particular interest, involved in the 
temporal evolution of a primordial CGW, in which the 
effective potential is constant for every t, namely 

R" SttG,^ , 

— = const = —Af2, (14) 

where M2 is a non- negative constant of units gr/cm^ [I]. 
In this case, Eq. (7) is written in the form 



+ k 



2 2 

' c — 



SttG 



0. 



(15) 



From Eq. (15), it becomes evident that, as long as 
M2 0, there exists a critical value of the comoving 
wave-number, 



SttG 



3c2 



■Mo 



(16) 



which discriminates the primordial GWs into oscillating 
{k > kc) and non-oscillating {k < kc) modes. According 
to [1] , such a property is met in a cosmological model with 
matter content in the form of a two-component fluid with 



Pir) 



M2 

i?2 



M4 
i?4 



(17) 



where M4 is an integration constant, also of units 
gr/cm?. 

For M4 = 0, Eq. (17) results in p R^'^ , which, by 
virtue of Eqs. (11) and (12), leads to p = —^pc^- This 
case, now known as the string- dominated Universe [19], 
was considered by Ford and Parker [15], in their effort 
for a natural renormalization of the energy-momentum 
tensor of the gravitational-perturbation field. 

A string-dominated Universe does not seem likely and, 
therefore, in this article, we consider M4 ^ 0. Accord- 
ingly, upon consideration of Eqs. (12) and (13), we dis- 
tinguish two cases, in Eq. (17), with respect to M2: (i) 
M2 = and (ii) M2 ^ 0. We consider each one of them 
separately. 



particle production [20]. Therefore, the choice of a solu- 
tion to Eq. (19) corresponds to the choice of an initial ( or 
final) quantum state for the gravitational-perturbation 
field. Eq. (19) admits the solution 



(r) 



irkc 



-ikcT 



(20) 



and the gravitational-perturbation field is in the well- 
defined adiabatic-vacuum state [21]. 



C. CGWs in the presence of cosmic strings 

Setting, in Eq. (17), A/2 7^ 0, corresponds to a radia- 
tion model contaminated by cosmic strings, in which, the 
two constituents (rclativistic particles and linear defects) 
do not interact with each other (cf. Eq. (13)). In this 
case, the constant effective potential (14), is related to 
the initial amount of the linear defects in the "mixture", 
A/2. Now, Eq. (15) results in 



(21) 



For every < fc < 00, the general solution to Eq. (21) 
can be expressed in the form 



TTjkc 



c\ e 



C2e 



(22) 



where ci and C2 are complex numbers, satisfying the 
Wronskian condition 



|ci| 



(23) 



(as a consequence of the mode normalization when the 
gravitational perturbations are quantized), and we have 
introduced the parameter 
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1/2 



(24) 



which measures the departure from the pure-radiation 
case, given by Eq. (20). 



B. CGWs in a radiation-dominated Universe 



III. A TALE OF COSMIC STRINGS 



For A/2 = 0, Eq. (17) reads 

A/4 



(18) 



In view of Eq. (12), this is the matter-energy content 
of a Universe filled with rclativistic particles - radiation. 
Now, the effective potential vanishes and Eq. (15) is 
written in the form 



k 



2„2, 



0, 



(19) 



i.e., it is a simple harmonic oscillator equation; like a 
wave propagating in flat space-time, where, there is no 



The presence of cosmic strings in the Universe is purely 
a question of topology [22]. In particular, after inflation 
(and reheating) the Universe enters in an early-radiation 
(ER) epoch [23], during which, the temperature drops 
monotonically (T ~ R~^). This cooling process may 
have resulted in the breaking of a fundamental U(l) local 
gauge symmetry, which, in turn, led to the formation 
of linear topological defects (for a detailed analysis see 
[7] and/or [S]). Within this era, the cosmological scale 
factor behaves as 



RERit) = ( 

Irr 



1/2 



Rer{t) = — , 

Trr 



(25) 
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where Tcr = 2 tcr marks the (conformal) time at which 
the cosmic strings are formed, and we have normalized 
R{tcr) to unity. 

By the time these hnear defects are formed, they are 
moving in a very dense environment and, hence, their 
motion is heavily damped, due to string-particle scatter- 
ing [24], [25], [20]. This friction becomes subdominant to 
expansion damping [24] at 

tcr ^ '^cr, (26) 

where ^ is the mass per unit-length of the linear defect. 
For r > , the motion of the long cosmic strings is essen- 
tially independent of anything else in the Universe and 
soon they acquire relativistic velocities [7] . Now, the gas 
of cosmic strings does not interact with the fluid of rel- 
ativistic particles and, hence, Eq. (17) holds. Therefore, 
we may consider that, for r > r*, the evolution of the 
curved space-time is driven by a two- component fluid, 
consisting of radiation (dominant) and cosmic strings 
(subdominant). Accordingly, marks the beginning of 
a radiation-plus-strings stage, during which, the Fried- 
mann equation (9) admits the solution [1] 

i?(r) ^ygsinh^M^,. (27) 

Nevertheless, the scale factor (27) can drive the Universe 
expansion only for a short period of time after , since, 
cosmic strings should (at any time) be a small proportion 
of the Universe matter-energy content. This means that 
the equation of state (17) has validity only for a limited 
time-period, otherwise, cosmic strings would, eventually, 
dominate on the Universe energy density (see, e.g., [19]). 

In fact, the radiation-plus-strings stage (if ever existed) 
does not last very long. Numerical simulations at the 
early 90's [27] - [30], as well as their present-time coun- 
terparts [.31] - [33], suggest that, after the friction be- 
comes unimportant, the production of loops smaller than 
the horizon, gradually, results in the scaling of the long- 
cosmic-string network. According to this process, the 
linear defects form a self-similar configuration, the den- 
sity of which, eventually, behaves as In this way, at 
some (physical) time tsc > t* the Universe re-enters in 
the (late) radiation (LR) era, 

RLBit) = Rsc ( — ) Rlr{t) = Rsc — , (28) 

before it can become string- dominated. The duration of 
the radiation-plus-strings stage is quite uncertain, mostly 
due to the fact that numerical simulations (which re- 
vealed the scaling of the long-strings network) can be 
run for relatively limited times. For example, the ear- 
liest treatments [29], [30] suggested that tsc — 30 t* 
{tsc — 5.5 r*), while, the most recent ones [32] raise this 
value to tsc — soot* {tsc — 17t»). According to Eq. (40) 
of [1], the former result leads to kccAr = kcc{Tsc — t,) = 
0.09, while, adoption of the latter result, would lead to 
fcccAr = 0.32. 



IV. GRAVITON CREATION 

Around r* and Tsc, both R{t) and R'(t) can be 
matched to be continuous (cf. [1]), but R"{t) is essen- 
tially discontinuous, acquiring a non-zero value through 
Eq. (14). Consequently, the scalar curvature of the 
space-time changes discontinuously both at r* and at Tsc- 
A discontinuous change in the scalar curvature produces 
gravitons [20], [34]. 

As a consequence, the evolution of the CGW modes 
in the transition of the Universe from an early-radiation 
epoch to the late-radiation era through a radiation-plus- 
strings stage, is modified as follows: 

• For r < the effective potential (14) vanishes. 
The temporal evolution of the metric perturbations 
is governed by Eq. (19), admitting the solution 
(20) and the gravitational-perturbation field is in 
an adiabatic-vacuum state. 

• For < T < Tsc the effective potential reduces 
to a non-vanishing constant (Ve// = k'^c^) and the 
evolution of the metric perturbations is driven by 
Eq. (21). Due to the discontinuous change of the 
scalar curvature at r*, during the radiation-plus- 
strings stage the gravitational-perturbation field 
is no longer in its vacuum state and the general 
solution to Eq. (21) is the linear superposition 
of positive- and negative- frequency solutions (22). 
Moreover, for every r, < t < Tsc, the primordial 
CGWs, being discriminated by the value of their 
comoving wave- number (in comparison to kc), do 
not evolve in a unique way: The GW modes with 
k > kc perform oscillations within the horizon, 
while, those of k < kc remain outside of the Hubble 
sphere and grow exponentially [1]. 

• Finally, for t > Tsc, Veff = and the temporal evo- 
lution of a metric perturbation in curved space-time 
is (once again) governed by Eq. (19). However, due 
to the discontinuity of the scalar curvature at Tsc, 
this time, the general solution to Eq. (19) is written 
in the form 



hr'ir) = ^-j-{ak e-^'^-- + e+'''^^) , (29) 

where the (Bogolubov) coefficients and /3fc should 
satisfy the Wronskian condition 

\a,.\^ ~ \Pk\' ^ I. (30) 

Eq. (29) determines the final quantum state of the 
gravitational-perturbation field. As a result of conformal 
invariance, this state is also an adiabatic vacuum, but, 
as long as /3fc ^ 0, it differs from the corresponding state 
of the early-radiation epoch. Clearly, the occurrence of a 
non-zero /?fc in the late-radiation era, would signal that 
/i™*(t) is not a pure positive- frequency solution, but con- 
tains also a negative- frequency component. This means 
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that, if the in-state of the gravitational-perturbation field 
is taken to be vaeuum, then, particles arc found in the 
out-state [21]. 

In what follows, we discuss the evolution of the modes 
(20), both with k > kc and with k < kc, in the transition 
of the Universe from the early-radiation epoch to the late- 
radiation era, through a radiation-plus-strings stage. The 
coefficients ci, C2, ak and 13k can be determined by the 
requirement that hkir) and its first derivative h'f.{T) are 
continuous across the boundaries and Tsc- 

Matching the modes (20) and (22), as well as their first 
derivatives, at t = t.^, we obtain 



and 



C2 = ^ } g-jfcc(7-Hl)r. 



2^/7 



(31) 



(32) 



for which, the Wronskian condition (23) holds. Similarly, 
matching the modes (22) and (29), as well as their first 
derivatives, at r = Tjc, we obtain 



and 



1 

47 



(7 + 1)^ g-'(7-l)fecAr _ _ 1)2 g+i(7-(-l)fccAr 

(33) 



1 



27 



sin(7A;cAr) e 



(34) 



where At = Tsc ~ and we have taken into account 
Eqs. (31) and (32). By virtue of Eq. (34), the number of 
gravitons [21] created in the mode denoted by < A: < 00, 
is written in the form 



|/3fc|2 = l, holds. 



Nk = l/^fel' = ^ (7' - 1)' sin2(7fccAr) (35) 

and the normalization condition |afe|^ 
The observable quantity Nk possesses a series of very in- 
teresting properties, depending on the value of k in com- 
parison to kc- We distinguish the following cases: 



A. Oscillating modes 



In this case, k > kc and Eq. (35) yields 



^4 sin ycAry^k'^ — k\ 
4fc2 P~~fc2 



(36) 



Obviously, for kc = 0, Nk>k^ vanishes, i.e., no gravitons 
are created in the absence of cosmic strings. 

On the other hand, for k ^ kc and/or At — >■ 0, Eq. 
(36) results in 



4A:2 



(37) 



Now, we distinguish the following cases: 



(i) At = 0: In this case, Nk>k^ = 0. Although this 
may look like an unexpected result, it is not. Usu- 
ally, the instantaneous transition from one epoch 
to another may lead to pathological results, espe- 
cially when high-frequency modes are concerned 
(see, e.g., [21]). However, in our case, as long as 
At = no transition ever takes place. The Uni- 
verse remains in the (early) radiation era, where no 
gravitons are produced. 

(ii) k = kc'. In this case, the number of gravitons cre- 
ated in the scaling of a long-cosmic-string network 
acquires the constant value 



Nk =^c'At\ 



(38) 



Eq. (38) suggests that, the number of gravitons 
in the ground state (kc) of the oscillating modes is 
well-defined and finite. 

(iii) Finally, for k 3> kc and k — > 00, Nk:>k^ ^ and 
particle production ceases. Consequently, there 
should be no ultraviolet divergences. 

On the other hand, according to Eq. (36), the expec- 
tation value of the number operator, Nk>k^, varies peri- 
odically with k. Hence: 

• For comoving wave-numbers satisfying the condi- 
tion 

c At^P - fc2 = mTT (to = 1, 2, ...), (39) 

^kyka ~ 0' admitting its absolute minimum value. 
The case m = is excluded due to the fact that 
Nk^ 7^ 0. As a consequence, creation of gravitons 
with comoving wave-numbers 

(to = 1,2,...), (40) 



(cAt)2 



never occurs by this mechanism, which (in 

quantum mechanics) is known as perfect transmi- 
tion (see, e.g., pp. 101 - 104). Clearly, the 
spectrum of these states is discrete and their po- 
sition on the fc-axis is determined (solely) by the 
choice of At 7^ 0. 

• On the other hand, the values of fc > kc that label 
the maxima of Nkyk^ obey the equation 



cAT^Jk^^ - kA cot (cAtV/c^ -kA = 2 



(41) 



which can be solved numerically, to give 

(n = 0,l,2,...). (42) 



49\^ 



36/ (cAt)2 

Once again, the position of these (discrete) states 
on the fc-axis is determined (solely) by the choice 
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of Ar ^ 0. Therefore, a possible detection of these 
waves, among other things, could reveal valuable 
information on the duration of the radiation-plus- 
strings stage, as well. 



dn^ = iJ^dijj j^-K^c" . Accordingly, the distribution of 
the gravitational-energy density, summed over the two 
polarization states of the gravitons present in the late- 
radiation epoch, is given by 



B. Non-oscillating modes 

According to [J], for < t < Tsc, the CGW modes 
with k < kc remain outside the horizon, where they 
grow exponentially. Nevertheless, from the quantum- 
mechanical point of view, low-frequency (non-oscillatory) 
modes are also created along the transition of the Uni- 
verse through the radiation-plus-strings stage, in the 
form of waves that tunnel through the barrier Veff of 
width At. Of course, doing so, they are subject to an 
additional exponential decrease (see, e.g., [21]), which, 
eventually, may counter-balance their original growth. 
Nevertheless, due to the tunneling effect, they can also 
contribute to the CGWs with k > kc created during the 
radiation-plus-strings stage, and, therefore, they should 
be taken into account in the late-radiation era. 

For k < kc, the expectation value (35) of the number- 
operator, Nk<k, = |/3fc<fcJ^, is given by 



<fcc 



1^4 sinh' 
4^2 



(cArv^ 



fc2 



k?. - B 



(43) 



Once again, either for fee = or for At = 0, Nk<cka 
vanishes, i.e., no tunneling ever takes place. On the other 
hand, for k — > kc, Eq. (43) is reduced to the constant 
value (38), verifying that, the number of gravitons in the 
state denoted by kc is well-defined and finite. 

In view of Eqs. (36), (38) and (43), it is now evident 
that, the number operator involved in the creation of 
gravitons along the transition of the Universe through a 
radiation-plus-strings stage is a continuous function of k, 
for every < fc < oo. 

Finally, for k <C kc, the expectation value of Nk<k^ 
reads 



3 

4fc2 



Nk<k^ - sinh^(fcecAT). 



(44) 



The apparent infrared divergence in Eq. (44) is fictitious. 
Indeed, these very low-frequency (very long-wavelength) 
waves are created way beyond the horizon [1]. As a con- 
sequence, they will enter inside the visible Universe at 
very late times, some of them (i.e., those with A: « 0) not 
even today. 



V. THE POWER SPECTRUM OF THE 
CREATED GRAVITONS 

During the (late) radiation stage of the cosmic ex- 
pansion, the physical frequency of a GW is defined as 
id — ck/R{T) and the number of states with frequen- 
cies in the interval to and uj + dtj is written in the form 



deg^ =r{uj)duj = 2hLu —^\pk\^ dio . (45) 
Once again, we distinguish the following cases: 

A. Oscillating modes 

In terms of fc > kc, Eq. (45) is written in the form 

d€k>k^ = V{k > kc)dk = (46) 
' (cAt^F"^^ 



1 he ,4 sin 



i?4(r) 47r2 



k'^-k?. 



k dk. 



Now, in view of Eq. (46), the corresponding power spec- 
trum [17] reads 



de 



V{k>kc) = ^^^^= (47) 

Bin^ (cAT^/k^ - kfj 



dk 

he 



4^2 Cp,V ; fc2_fc2 

where we have introduced the (physical) critical wave- 
number at every t, as 



kc,Ar) = kc/R{T). 



(48) 



For k close to kc, V{k > kc) ^ k, while, in the large- 
limit, it results in a damped oscillation 



r(k> kc) sin2(A:cAT). 
k 



(49) 



According to Eq. (49), the production of gravitons with 
k > kc is suppressed at high frequencies. 

Similarly, the logarithmic spectrum is written in the 
form 



de 



d(lnfc) 



kVik > kc) 



(50) 



in^ ( cArJk^ - k^ 



he sill 
47r2 



A:2 - fc2 



Eq. (50) becomes particularly transparent (and useful) 
in terms of the dimensionless parameter (see, e.g., [1]) 



kc 



(51) 



which measures the comoving wave-number of the oscil- 
lating modes in units of their ground-state counterpart, 
kc. In terms of x, Eq. (50) is written in the form 



de 



k>ka 



d(lnk) 



kV{k > kc) = xV{x > 1) 



(52) 



he 4^ sin^ [kccAT\/x^ - 1) ^ 
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or, else, 

xV{x>l) ^ ATT^ikccArf 2hckc^^{T) ^ 



sin^ [kccAr^/x'^ - l) 
(fcccArVa;2 - l)^ 



(53) 



where Acp,^ = 2Tr /kcp,^ is the correlation (wave)length, 
associated to kc„u . 



Cl 




100 150 

(k/kc) 



250 



FIG. 1: The logarithmic spectrum k Vik), of the high- 
frequency part (k > kc) of the metric perturbations created in 
the scaling of a long-cosmic-string network (normalized over 
47r tc), versus k in units of kc, for Ar = 5.5 r» (according 
to [1]). 

At a fixed t > Tsc, the gravitational energy contained 
within the volume Ai? , is 2hckr ^, i.e., the energy of a 
graviton with physical wave-number kc ^ summed over 
the two polarization states. Therefore, the quantity 



ec(r) = A^fe, 2?icfccp,(T) 



(54) 



(where A'^^^ is given by Eq. (38)) represents the energy 
density of the gravitons created in the (ground) state de- 
noted by kc, at every r > Tsc and the logarithmic spec- 
trum (53) reads 



xr{x > 1) = 47r edr) 



sin^ (fcccAr-v/ x'^ — l) 



(55) 



We can easily verify that, for a; = 1, one is left with the 
energy density of the gravitons in the state k = kc (within 
the solid angle of a sphere - 47r), while, as we depart from 
the lowest-allowed comoving wave- number of the oscillat- 
ing modes, the gravitational power is distributed among 
the various fc-intcrvals according to Eq. (55). In partic- 
ular, on the approach to kc, the logarithmic spectrum is 



quadratic in k, while, for k ^ kc {x ^ oo), it reduces to 
a periodic function of the comoving wave-number with 
constant amplitude (Fig. 1) 



dek>k^ 
d{\nk) 



sivl' [kcAr) . 



(56) 



The importance of this result rests in the fact that the 
logarithmic spectrum is a re-scaling of the function 



1 dtg.^ _ ^ _ 1 dt 
Ecr d{h\k) Ecr d{\n J) 



gw 



(57) 



(where Ecr is the critical- energy density for closing the 
Universe), which characterizes the intensity of a stochas- 
tic background oi CGWs [■-)(»]. 

Therefore, if a relic CGW signal has been filtered by 
a long-cosmic-string network [1], its profile would have 
been modified along the lines of Eq. (55) (or Eq. (56)), 
thus, resulting in a periodic function of the frequency, 
the period of which depends solely on the duration of the 
corresponding radiation-plus-strings stage (Fig. 1). 



B. Non-oscillating modes 



In this case, in terms of the dimensionless parameter 
X ~ 7^ (< 1), the distribution of the gravitational-energy 
density of the waves with k < kc, which are present in 
the late-radiation epoch, is given by 



rfefe<fe, = 47r edr) 



sinh" (/cccArVl — a;-^) 
(fcccArVl - a;2)^ 



X dx, (58) 



where eclr) is given by Eq. (54). Accordingly, the corre- 
sponding logarithmic spectrum reads (notice that, here, 
Vk stands for of [1]) 



dek 



<fce 



d{\nk) 



kVik < kc 



(59) 



47r edr) 



sinh (fcccAryl 



(kccArVT 



and is presented in Fig. 2, which, in fact, corresponds to 
a (magnified) "detail" of the low- left corner of Fig. 1. 

Accordingly, if the Universal evolution has undergone 
through a radiation-plus-strings stage, then, in connec- 
tion to [I], the anticipated profile of a stochastic GW 
signal should be modified as follows: 

• Every GW background of cosmological origin that 
created earlier than (e.g., the metric perturba- 
tions predicted by inflation), is being filtered by the 
long-cosmic-string network, thus, being discrimi- 
nated to (high-frequency) oscillating modes and 
(low-frequency) non-oscillating ones. The separa- 

trix between these modes is set at kc = ^^^'h, 
where M2 is the initial amount of the linear de- 
fects, corresponding to a present-time frequency, 
fc, given by Eq. (56) of [1]. 
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Q. 




parameter). According to [1], we have 

-2 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

(k/kc) 



FIG. 2: The logarithmic spectrum kP(k), of the CGW modes 
with k < kc, i.e., those created outside the horizon (normal- 
ized over 47r Ec), versus k in units of kc, for At = 5.5 r* (solid 
line) and At = 17 r, (dashed line). 



• Due to quantum-gravitational effects, the high- 
frequency part [k > fee) of the corresponding CGW 
signal is being further amplified by a factor of 100 
(see, e.g., Fig. 1), while, its profile is turned to 
a periodic function of the frequency, the period of 
which depends (solely) on the duration (Ar) of the 
radiation-plus-strings stage. To the best of our 
knowledge, such an amplification mechanism of a 
stochastic GW signal has not yet been considered 
elsewhere. 



• According to our model, the only interference 
between an early-radiation epoch and the late- 
radiation era comes from a period of significant 
presence of cosmic strings (encapsulated in the cos- 
mological model). Therefore, it is natural to as- 
sume that any energization of the gravitational- 
perturbation field (graviton production) is due 
to a corresponding energy-loss of the long-linear- 
defects network. In this context, it is possible 
that the actual scaling process might have been 
completed a little bit earlier than what is cur- 
rently anticipated by the scientific community (see, 
e.g., [29], [30], [32]). 



tc > 



C2 J 



tr 



(60) 



where tcr is the time at which the Universe acquires the 
critical temperature below which the cosmic strings are 
formed (for GUT-scale strings, tcr ^ 10~^^ sec) and 7* 
is the number of correlation lengths (4 < 7* < 7) inside 
the horizon at t^ = (^)~^tcr (see, e.g., [7]). 

The dimensionless quantity [Gfi/c^) characterizes the 
strength of the gravitational interactions of strings. The 
observations give an upper limit on the value of this pa- 
rameter. In particular, the current CMB bound is [37] 
(Gfi/c^) < 1.3 X 10~^ while, more recent studies on a 
cosmic-string contribution in the WMAP data [38] have 
yielded the tighter bound (G/i/c^) < 3.3 x lO"'^, some- 
thing that is also confirmed by gravitational lensing ob- 
servations: {Gfi/c^) > 4 X 10"'^ [39], [40]. Therefore, in 
what follows, we may consider (Gn/c^) = ^ x 10~^, so 

that, tc > 10^2 tcr ^ 10"^^ sec. 

As far as the oscillating modes are concerned, kc is the 
lowest-allowed comoving wave-number. Therefore, after 
tc, no oscillating CGW modes (created by this mecha- 
nism) are left, to enter inside the Hubble sphere. For 
t > tc, the only contribution to the CGW-energy density 
comes from the low-frequency, non-oscillatory modes, of 
comoving wave- number k < kc- 



VI. THE EXTRA GW-ENERGY DENSITY 

Now, we may proceed even further, to calculate the 
(integrated) extra GW-energy density present in the late- 
radiation epoch, due to the metric perturbations created 
in the scaling of the long-cosmic-string network. Accord- 
ing to what previously stated, two kinds of modes are 
present in the late-radiation era. We consider each case 
separately: 



A. Oscillating modes 

The integrated energy density of the oscillating modes 
{k > kc), created in the transition of the Universe 
through the radiation-plus-strings stage, is given by 



^gw {k > kc) 



dek>kAk) 



dck>kA^), (61) 



As we have already mentioned, the number of gravi- 
tons, Nk^ , created in the ground state of the oscillating 
metric perturbations, is well-defined and finite. These 
modes enter inside the horizon at some (coordinate) time 
tc > tsci at which, their physical wavelength (defined as 
\h{tc) = ^cR{tc)) becomes smaller than the correspond- 
ing Hubble radius (initc) = c/H{tc), H being the Hubble 



which, in view of Eqs. (38) and (55), is written in the 
form 



£gw{k > kc) = TT ec(T) 



1 

°° sin^ {2./W,^/x^^ 



(62) 



xdx 
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or, else, 



1 



£gwik > kc) = 47r ec(T) — — x 



°° sin^ (y]Vfc7Vx2 - 1) cos2 (yiVfe7Vx2 - 1) 



1 



(63) 
xdx. 



In fact, the high-frequency gravitons, responsible for the 
infinite upper- limit, can be produced only by an instan- 
taneous change in the scalar curvature. In our case, 
the highest frequency (which can be produced) is deter- 
mined by the speed of the transition (the duration of 
the radiation-plus-strings stage), otherwise particle cre- 
ation is adiabatically suppressed (see, e.g., [20], [21], [34]). 
Hence, 



jkc = c\l k"^ — fc2 < — — 
At 

(cf. Eqs. (21) and (22)) and, therefore. 



(64) 



(65) 



Now, in terms of the parameter x, Eq. (63) is written in 
the form 



£gw{k > kc) = 47r ec(T) 



(66) 



X 



dx 



and the infinite upper-limit can be recovered under the 
assumption N^^ ^ 1 (in fact, this assumption is quite ac- 
curate, since, according to [1], Nk^ = 2x 10^'^). However, 
in this case, we have 

sin(y]V^x) = <X-^</V + 0(</^) 
cos (v/iVfeTx) = 1 - ^iVfc^x' + O (Nl) . (67) 

Accordingly, discarding terms of order O (7V^ ) , Eq. (66) 
is written in the form 



£gw{k > kc) ~ 47r edr) 



thus, resulting in 



X-^Nk,x^ ] dx, 



(68) 



5tt 



1 



£gw{k > kc) — —— Cc{t) — 



12 ' Nk^ 

By virtue of Eq. (54), Eq. (69) takes on its final form 

1 5hc 



(69) 



^gwik ^ kc) 



kt 



R'^{t) 487r2 " 



(70) 



Inside the Hubble sphere of the late-radiation era, the 
energy density of the gravitons with k > kc which were 



created at earlier times tsc < t' < t and the wavelength of 
which became less than the horizon due to the expansion 
of the Universe at tsc < t < tc, is given by the sum of all 
the previous increments 5£gw{t') being redshifted by the 
factor R'^{t')/R^{t) 

eg.Ats.<t<tc) = l^^[-^^.^)'^t'^in) 

egAts.<t<tc) ^ ^<.(u(^)'ln(^ 

For t > tc (i.e., deep inside the late-radiation era), this 
constant amount of energy density is being (further) red- 
shifted due to the cosmological expansion. Therefore, 



(t) = £gw {tsc <t <tc) 



R\tc) 

RHt) 



.k>k^/,\ _ 5?IC 4 ftc\ , t tr 



247r2 



while, for t > tc, there are only the non-oscillatory modes 
available to enter inside the horizon and contribute to the 
extra GW-encrgy density. 



B. Non-oscillating modes 

The extra GW-energy density (present in the late- 
radiation era) of the modes with k < kc, created in 
the transition of the Universe through the radiation-plus- 
strings stage, is given by 

Egw{k<kc)= I d€k<kAk) ^ / dekKkA^)-, (73) 
Jo Jo 

which, by virtue of Eqs. (38) and (58), is written in the 
form 



£gu,{k<kc) = tt€c{t)-^ 



Nk^ 

1 sinh^ (2y]VA~Vl - 



(74) 



1-X2 



' dx 



Now, in terms of the parameter = \/l — x^, Eq. (74) 
reads 



£gw{k < kc) = An e^r) 



(75) 



1 sinh^ (V^0) cosh' (Ta^ 



which, in the limit Nk^ ^ 1, results in 



£g^{k < kc) ~ An e,(r) / U + -NkA' d<P, (76) 



yielding 



£gw{k < kc) ~ 27r edr) 1 + -Nk 



(77) 
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or, else, 



{k < kc) 



1 he 4 



(78) 



For the same reasons that led to Eq. (71), the energy 
density of the gravitons with k < k^ which were created 
at earlier times tc < t' < t and the wavelength of which 
became less than the horizon due to the expansion of the 
Universe at tc < t < tree (where tree is the recombination 
time), is given by 



he 



e^<^^ii)^N,-kl^{te) 



In 



(79) 



C. The overall energy density 

With the aid of Eqs. (72) and (79), we can deter- 
mine the (overall) extra GW-energy density, due to the 
metric perturbations created in the scaling of a long- 
cosmic-string network and entered inside the horizon at 
tc < t < tree- Wc obtain 



(i > tc) 



(80) 



oscillating 



non — oscillating _ 



where we have marked the contribution of both the os- 
cillating and the non-oscillating modes. By definition, 
the physical wavelength of the mode kc at tc equals to 
the horizon length (at tc this mode enters into the visible 
Universe), so that 



kc„h (tc) 



2-K 



27r 2TiH{tc) 



Kuitc) initc) c 
Accordingly, Eq. (80) is written in the form 



(i > tc) 



24 V tsc 



(82) 



^fc. In I - 



while to — 4.3 x 10"'^'' sec denotes the present epoch. 
According to Allen [20], during the inflationary epoch, 
Hinf = 8 X 10'^'* see~^, while, within the late-radiation 
era, Her = (2tcj)~^ = 5 x 10'^° sec~^. Finally, admitting 
that At = 5.5r* (see, e.g., [1]), we obtain Nk^ = 2x 10~^, 
tsc — 30 and, therefore, tc — 10^ tsc- Inserting these 
values into Eq. (83) wc obtain 



^rps 



^it>tc) 



4.28 X 10^ 



-gw 



secy 

t) 



1 + 10"^ In 



t 



(84) 

We observe that, at late times, the only contribution to 
comes from the second term in the brackets on the 
right-hand- side of Eq. (84), i.e., from the (low-frequency) 
non-oscillatory CGW modes (a not unexpected result). 
Nevertheless, within the late-radiation era, this term in- 
troduces only minor corrections to the p- behavior of the 
CGW-energy density. 

It is generally admitted that the main source of relic 
gravitational radiation is the inflationary amplification of 
metric perturbations, taken place as early as 10""^'* see 
after the Big Bang (see, e.g., [12]). Quite later, t >tsc ''-^ 
10~^^ sec, after being filtered by a long-cosmic-string 
network, the original number of the relic (inflationary) 
GWs is further increased by an additional amount of 
CGWs (cf. Eqs. (36) and (43)), due to the quantum- 
gravitational creation of gravitons along the transition of 
the Universe through the radiation-plus-strings stage. In 
this case, Eq. (84) indicates something very interesting: 
At times tc ~ 10~^^ see < t < 6540 sec, the contri- 
bution of these gravitons, e^^", to the (already) existing 
gravitational-energy density , is (much) larger than 
efj itself. 

Certainly, at late times, the amount, Sg^, of the ex- 
tra gravitational-energy density (created in the scaling 
of a long-cosmic-string network) becomes negligible, as 
a consequence of the cosmological expansion. Neverthe- 
less, it may have affected significantly the cosmological 
evolution within the first (and most interesting) 109 min 
of the history of the Universe, thus leaving imprints on 
the anticipated stochastic GW power spectrum. Accord- 
ingly, for f > fc (see, e.g., [1]), the spectral intensity fig^, 
of the CGW signal is turned to a periodic function of the 
frequency, while being amplified by (almost) two orders 
of magnitude (Fig. 1). 



which, in view of the Appendix A, results in 



^rps 



^{t>tc) 



247808 



81 



2 4 
-TT 7* 



n1 





Her 


2 






. Hinf . 







24 



■In 



Nk^ In 



(83) 



where e^J^ is the GW-energy density predicted by (de 
Sitter-type) infiation, and Her are the values of the 
Hubble parameter at infiation and at tcr, respectively. 



VII. CONCLUSIONS 

Assuming that there was a period in the early Universe 
in which the matter content could be modelled by a two- 
component fiuid (consisting of radiation and a cosmic- 
string fiuid), we demonstrate that the existence of such 
a (radiation-plus-strings) stage would have resulted in 
a measurable effect on the CGWs' power spectrum. In 
fact, along the transition from an early-radiation epoch to 
the late-radiation era through the radiation-plus-strings 
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stage, gravitons with a highly-rccognizablc profile eould 
have been produced. 

As a consequence, at high frequencies, the CGW spec- 
trum would no longer be scale invariant, but, rather, a 
•periodic function of the frequency, the period of which 
depends (solely) on the duration of the radiation-plus- 
strings stage. 

But, what is more important for the detection of 
CGWs, is that, the extra amount of metric perturba- 
tions created in the scaling of a long-cosmic-string net- 
work would lead to the amplification of the relic GW 
signal by (almost) two orders of magnitude. 

The extra- energy density of the CGWs created in the 
transition from an early-radiation epoch to the late- 
radiation era through the radiation-plus-strings stage, 
not only is added to the corresponding inflationary quan- 
tity. In fact, at times 10^^^ sec < t < 6540 sec, i.e., 
within the first (and, certainly, most interesting) 109 min 
of the Universal evolution, dominates over every other 
form of GW-energy density. 

According to our model, the only interference between 
an early-radiation epoch and the late-radiation era comes 
from a period of a significant presence of cosmic strings 
(encapsulated in the cosmological model). Therefore, 
it is natural to assume that any energization of the 
gravitational-perturbation field (graviton production) is 
due to a corresponding energy-loss of the long-linear- 
defects network. In this article we have neglected the 
backrcaction of the created gravitons on the curved cos- 
mological background (it will be the scope of a future 
work); nevertheless, we expect that the actual scaling 
process might have been completed a little bit earlier than 
what is currently anticipated by the scientific community. 
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Appendix A 

The contribution of the gravitons created in the tran- 
sition from an early-radiation epoch to the late-radiation 
era through a radiation-plus-strings stage to the total 
amount of the gravitational-energy density in the Uni- 
verse, is given by 

sll^it>t^)^16n'^H\t^) (I'' 



(Al) 



(cf. Eq. (82)), which, upon consideration of the Fried- 
mann equation 



reads 



ff'(ic) = ^e™d(ic), 



e7:{t>t.)^^^il\t,)ep,e,^a{t) 



24 ^""(i, J ^^'^'^ ^""(c 



(A3) 



where tpi = 5.44 x lO""*^ sec is the Planck time. Taking 
into account that, in terms of Sraditc)-, the mass-energy 
density in the matter-dominated era is given by 



£matt{t) — £rad{tc) 



where to denotes the present epoch, we find that 
e7J{t > t^) = ^iJ^(tc) t% e^Mt) 



(A5) 



which, upon definition of as the Hubble parame- 

ter at (dc Sitter) inflation, can be written in the more 
convenient form 



2048 



Hit,) 



inf 



£matt{t) 



to 



24 ^""(t, J ^^^^^''(i. 



(A6) 



and, by virtue of Eq. (4.10) of [20] (or Eq. (16) of [::i4]), 
takes on its final form 



2048 4 

7r 

g-^ 9 



Hit,) 



H 



inf 



iA7) 



Eventually, taking into account that, in a FRW radiation 
^, with the aid of Eq. (60), Eq. 



model we have 



H(ta 



(A7) resuhs in Eq. (83). 
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